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1 Introduction. 



Our research was targeted to analyze, through computer simulations, the res- 
olution power of a spectrometer with different kinds of magnetic fields and at 
different energies. 

Throughout our analysis, we have focused our attention on the tt~ tagging in 
the A — > pvT reaction. 

We wanted to study the possibility of making a good tagging at high energies 
(> fOOGey), i.e. a good A momentum reconstruction with the software at our 
disposition. 

To this end the work has been divided in two parts: 

• The experimental setup definition. 

• The proper simulation. 

2 Definition of the setup. 

In line of principle the most straightforward definition of an experimental setup 
of this kind must concentrate on three points: 

1. The target. 

2. The magnet (i.e. the magnetic field). 

3. The detectors. 

2.1 The target. 

At our level of analysis the presence of a real, material target was deeemed 
unnecessary, so that, to simplify the study and not to introduce in it too many 
variables, we have decided to identify it as a geometric point in a first time and 
with an area later. 

2.2 The magnet and the field. 

To greatly simplify the first part of the analysis and to be able to use quick and 
simple analytic formulas for the motion of particles in a magnetic field, we have 
decided to define initially the field as homogeneous of dipolar type with: 

Bx = ^ By = ^ Bz = const 

The constant value of the field has been set to 1.5 or 3 T, according to the 
momentum of the A. 
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The study of the effect of more complex and more interesting field configu- 
ration has been delayed to a later time, when the problem had been correctly 
mastered in its easier configuration. 

2.3 The detectors. 

The detectors' dimensionment and placing was certainly the more complicated 
task of the first part of this study, because their optimization depends closely on 
the chosen magnetic field, and we wanted to modify this field during the study. 
At the same time there is also an energy dependence of the positionment and of 
the resolution: at different value of the A momentum (pa), different resolution 
are required to obtain the same precision in the momentum reconstruction and 
different detector positions are required to take in account the different mean life 
of the A. 

To this end and because of the necessities of our momentum reconstruction 
program, we have decided to put an high number of detector of the same stan- 
dardized form. This form is next to the optimum only for the dipolar field, but 
we think that the generality of our analysis can go along with the lack of some 
optimization. 

The positionment of the detectors was made considering the particle motion 
in the uniform field at different momenta in a range from 0.5 to more than 100 
GeV/c. 

Such a wide range forced us to a compromise among the best possibilities 
present at lower (< lOGev/c) and higher momenta. The final setup is intended 
to perform decently for momenta up to 250 GeV/c, but clearly is not the most 
efficient or cheaper that can be envisaged for a particular momentum subrange 
or a different kind of field. 

This setup is showed in 1 within the reference system used for the simulations. 

The detectors 1-3 have a resolution of 250/im, the detectors 5 - 8 of 0.5cm 
and the others of 1cm. 

The detector distances from the target (origin) and their angular openings 
are summarized in Table 1. 
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Figure 1: Detector position in the reference system of the simulations. The 
dummy target is the small square. The object before the target is a catcher 
intended to stop low momentum particles from spiralizing around the target. 
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TABLE 1 - Detectors positions and openings. 



DETECTOR NUMBER 


RADIUS (cm) 


(f) (degrees) 


1 


10 


2 * 28 


2 


30 


2 * 28 


3 


50 


2 * 28 


4 


100 


2 * 27 


5 


120 


2 * 27 


6 


140 


2 * 27 


7 


160 


2 * 27 


8 


180 


2 * 27 


9 


210 


2 * 26.5 


10 


230 


2 * 26.5 


11 


250 


2 * 26 


12 


260 


2 * 26 


13 


270 


2 * 26 


14 


280 


2 * 26 


15 


290 


2 * 26 



As illustrated in 1, a 0.5 - 1 degree open window has been left in the more 
distant detectors to allow very high momenta to pass through without adding to 
the background. 

Eor the purposes of definition in the simulation software, the thickness of the 
detectors has been specified in 1 mm, but this hasn' t a big relief, because the 
detectors are considered as if they were made of vacuum, that is not interacting. 
This simplification was done in order to be able to evaluate the effects of the 
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different fieids without superimposing a random scattering component. 



3 The simulation. 

For the computer simuiation we have used mainiy two programs: a Monte Cario 
simuiation program based on the GEANT paclcage [f ] and a momentum recon- 
struction program essentiaiiy based on the MOMENTM routine [3], which was 
inciuded up to some time ago in the CERN iibraries. 

The physicai probiem on which we focused our efforts was the determination 
of the momentum of a decaying A particle through the detection of the generated 
pion alone. 

This technique (" tagging ") has the advantage of allowing one to deal with 
lower momenta, in fact the pion carries away at most approximately one fourth 
of the momentum of the decaying A. So the momentum reconstruction becomes 
easier because the magnetic rigidity is reduced in the same ratio and the angle 
of emission at the decay increases. 



3.1 Kinematics. 

In a two body reaction we can write [5], in the lab reference system: 



2 2 2 2 / ^ \ ^ ^ 

2pip3Cos63 = — — m2 — + 2[Ei + m2)E3 — 2Eim2 

in our case (A decay) we have that subscript 1 refers to the A, particle 2 is 
non existent, 3 is referred to the pion and 4 to the proton, so that we can simplify 
to obtain: 

2pAP7rCos9x7r = - - + 2EaE^ 
if we remember that : 

/I PiTX 

COS Att = 

Ptt 

if we assume that p^ is parallel to the x axis, with a bit of algebra we can 
obtain for the A momentum the following expression: 



-Vp^, ± E^^V^ + mi(p2^ - El 
= E^-p . 



where 



As the ^Att angle becomes smaller with the growth of pA, the reconstruction 
of the X component becomes more and more critical : a small error can lead 
the square root in equation 1 to become imaginary. 

Other formulas have been tried, but they behave in no better, or even worse 
manner. 

3.2 Structure of the simulation. 

The Monte Carlo simulation program is based on GEANT, it generates a A 
particle in a predetermined zone of space, generates the decay, then tracks it 
according to the laws of kinematics and of motion of a charged particle in a given 
magnetic field. 

Hits of the particles with the simulated detectors are registered. This recorded 
hits are feeded as input in the momentum reconstruction program, which outputs 
the momentum module and components of the concerned particle. This particle, 
in our tagging, is the pion, knowledge of these data allows us to reconstruct 
through the 1. 

The MOMENTM routine, which is the core of our momentum recostruction 
program, gives its output (momenta) in the position corresponding to the first 
space point in input. We of course do not know the decay vertex position, because 
it is ruled by randomness through the low of particle decay, but we must have 
the X component of the pion momentum in the vertex with the highest possible 
precision, because the formula 1 is very sensitive to this valule, especially as the 
energy grows and the decay angle decreases. So we have made a modification 
in the original MOMENTM routine in order to be able to find the decay point, 
unfortunately the method employed is rigorously valid only if the A momentum 
vector is perfectly coincident with one of the axes of our reference system. This 
condition is easily satisfied if the A is generated in the origin with the required 
direction of the momentum vector, but otherwise it is violated, so we will able 
to reconstruct with the best accuracy only in this particular case. If the quoted 
condition is not satisfied, the vertex finding method still gives meaningful ap- 
proximate results, but they deteriorate quickly the more we are far from the best 
situation, and above all with the rising in energy. 



We have done our simulation with three kind of magnetic field, whose analytic 
expressions are : 



4 The results. 



B 



X 



, B, 



y ~ 



^ Bz = const 




for the homogeneous field. 
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X 

Bx = ^ By = ^ Bz = const x atan[—) (3) 

y 



and 



where 



„ „ „ const , , 

' ^ r + 20 ^ ' 



r = 

10 

if the coordinates are in centimeters and the field in kilogauss (GEANT units). 

We made two main sets of simulations, one with Pa = 5 GeV/c and the other 
with Pa = 250 GeV/c^ the constants were calculated to give B^ax = 1-5T in 
the first case and B^ax = S.OT in the second. 

In both cases we made separate simulations for A generated in and outside 
the origin of the reference system, always with parallel to the X 9jX1S. 



4.1 Runs at 5 GeV/c. 

With the field 2 we obtained 0.1 % both for particles starting in the 

origin or a few centimeters outside. 

With the field 3 we had 0.3 % again for both kinds of A generation. 

With the field 4 the results 10 % for a A starting 1.5cm out of 

PA ° 

the origin and ~ 3 % for the well behaved case. 

° PA 

4.2 Runs at 250 GeV/c 

Equation 1 is now much more critical and we had 3.5 % for the field 2 

... 
and the A generated in the origin. 

If the the generation occurs at r = 0.15cm we still have sa 4 %, but if 

. . . PA 

r = 0.3cm, the square root in 1 becomes immaginary and our formula fails. 

With the fields 3 and 4 we had ^ sa 35 % and ^ sa 75 % respectively, 
already in the case of the A starting from the origin. 



5 Conclusion. 

Our simulations have shown that the tagging is a usable and useful tecnique at 
momentum range of a few GeV/c even in a inhomogeneous field. On the other 
side it is very difficult to apply beyond a few tens of GeV/c. 

Actually we made a few simulations for = 20G'ey/c and for p^ = lOOGeV/ c. 
In the first case the resolution on p^ worsened of a factor of at least two, with 
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respect to the bGeVj c runs, while in the second case the situation was almost 
identical to that of the 2bQGeV j c runs. 

So we can conclude that without a more refined and specialized software 
for the momentum reconstruction (and perhaps refined and specialized detector 
setup) the tecnique illustrated for this study appears better suited for low and 
intermediate energies experiments. 
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